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We investigate pliantom models with power-law potentials and we extract the early-time, 
"tracker", solutions under the assumption of matter domination. Contrary to quintessence case, 
we find that energy positivity requires normal power-law potentials instead of inverse power-law 
ones, with the potential exponent being bounded by the quadratic form. In addition, we analyti- 
cally present the general cosmological solution at intermediate times, that is at low redshifts, which 
is the period of the transition from matter to dark-energy domination. The comparison with the 
exact evolution, arising from numerical elaboration, shows that the tracker solution agrees with the 
later within 2% for redshifts z > 1.5, while the intermediate solution is accurate within 2% up to 
z r; 0.5. 
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I. INTRODUCTION 

Recent cosmological observations support that the uni- 
verse is experiencing an accelerated expansion, and that 
the transition to the accelerated phase has been realized 
in the recent cosmological past |l|. In order to explain 
this remarkable behavior, and despite the intuition that 
this can be achieved only through a fundamental theory 
of nature, physicists can still propose some paradigms for 
its description, such are theories of modified gravity , 
or "field" models of dark energy. The field models that 
have been discussed widely in the literature consider a 
canonical scalar field (quintessence) [1, a phantom 
field, that is a scalar field with a negative sign of the 
kinetic term [sf, or the combination of quintessence and 
phantom in a unified model named quintom Q. 

In field dark energy models one can find potential- 
independent solutions [7| , but in general the cosmological 
evolution depends significantly on the potential choice. 
In the quintessence case, one well studied potential is the 
inverse-power law one [1, H, Q . This potential exhibits 
"tracker" behavior, that is for a large class of initial con- 
ditions the cosmological evolution converges to a common 
solution at late times [l^]- Furthermore, in such models 
the quintessence energy density remains small at early 
and intermediate times and thus the known cosmological 
epochs are not affected. Finally, a theoretical justifica- 
tion of power-law potentials can arise through supersym- 
metric considerations [11.]. 

The cosmological evolution in such quintessence mod- 
els at early times, is dominated by the matter-fluid and 
one can obtain the tracker solutions analytically [l^ . At 
intermediate times the quintessence field is still small, 
but non-negligible, and one has to rely on perturbative 
analytical expressions 0]. However, such an extension 
to low redshift is very useful, since observations like su- 
pernovae la, WMAP and SDSS ones, are related to this 
period [ij. 



In this work we are interested in investigating the 
phantom scenario in power-law potentials, both at high 
and low redshifts, and provide the tracker solutions and 
the perturbative analytical expressions respectively. Due 
to the "inverse" kinetic behavior of phantoms in potential 
slopes, we expect that we have to consider normal power- 
law potentials instead of inverse power-law ones. Indeed, 
starting with an arbitrary power law, we find that phys- 
ically meaningful cosmological evolution excludes nega- 
tive exponents. 

The plan of the work is as follows: In section |TT] we con- 
struct phantom models in power-law potentials and we 
present the matter-dominated solutions. In section lTTTl we 
derive the cosmological solutions at intermediate times, 
that is when the phantom dark energy is non-negligible, 
but still sub-dominant comparing to the matter content 
of the universe. In section IIVI we compare our analyti- 
cal expressions with the exact numerical evolution of the 
phantom scenario. Finally, in section |V] we provide some 
applications of the model discussing its cosmological con- 
sequences, and we summarize our results. 



II. COSMOLOGICAL EVOLUTION DURING 
MATTER DOMINATION 

Let us construct the simple phantom cosmological sce- 
nario in a flat universe. The action of a universe consti- 
tuted of a phantom field (/> is @ : 



5 = 



4 + V{cp)+C, 



,(1) 



where we have set SttG — 1. V^(0) is the phantom field 
potential and the term accounts for the (dark) matter 
content of the universe, considered as dust. Finally, al- 
though we could straightforwardly include baryonic mat- 
ter and radiation in the model, for simplicity reasons we 
neglect them since we are interested m. z < 2Q era. The 
Friedmann equations are 
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(3) 



while the evolution equation for the phantom field is: 

p^ + 3Hip^+p^) = 0, (4) 

where H = a/a is the Hubble parameter and a — a{t) the 
scale factor. In these expressions, pm is the dark matter 
density, and p^ and Prj, are respectively the density and 
pressure of the phantom field given by: 



(5) 
(6) 



Therefore, we can equivalently write the evolution equa- 
tion in field terms as: 



(7) 



The system of equations closes by considering the evo- 
lution of the matter density, which in the case of dust 
reads simply pm = Pmo/O'^, with p„io its value at present. 
Finally, the dark energy equation-of-state parameter is 
given by 



W : 



(8) 



As mentioned in the introduction, in the present work 
we study phantom evolution in power- law potentials. 
That is we consider: 



(9) 



where k is a constant with units of ra^^" . Note that 
we prefer to insert the minus sign in the exponent in 
order to coincide with the usual notation of the literature. 
A priori we do not restrict the value of a. However, 
as we are going to see, energy positivity will force a to 
be negative and bounded, that is potential ([9]) will be a 
normal (and not inverse) power law. 

It proves convenient to express the aforementioned cos- 
mological system using the scale factor a as the indepen- 
dent variable, since it is straightforwardly related to the 
redshift z which is used in observations. Following Q we 
define 



eV da) 



^Pm + P4,) 6V da. 

and thus the first Friedmann equation gives 

{Pm + V), 



l-x 



(10) 



(11) 



where V stands for V{(j)). Thus, one can simply write: 

xp„i + V 



1 



P4> 



XPm ~ V{1 - 2x) 



Xpri 



+ v 



- X 

V{1 - 2x) 



Xpr, 



V 



(12) 
(13) 
(14) 
(15) 



Comparing to the quintessence case, and apart from 
the obvious expectation that x is now negative, we ob- 
serve, that it is also bounded. In particular, requir- 
ing phantom total energy positivity (p^ > 0) we obtain 
X > -V/i'iH'^). Finally, using expressions ([T1|)-(I14]), the 
field evolution equation ([7]) becomes 



2 J," 



acj)' 



-(5-3.)-.^^-^)^^^'^ 

2 Pr, 



■V 



'dSJ 



0, (16) 



with prime denoting the derivative with respect to a. 
This equation is exact and accounts for the complete dy- 
namics of the phantom scenario. 

In this section we desire to present the tracker solu- 
tions during the matter-dominated era. Thus, we con- 
sider Ip^I ^ or equivalently |a;| <C 1. In this case 
equation (|16p is simplified as: 



2 J," 

a 0, 



Sac 



(0) 



■(0) 



Pm acj) 



(17) 



= 0(0) 



The zero subscript in parentheses 



with ^ = % 

a<p aq> 

denotes just this zeroth-order solution in terms of x, or 
equivalently in terms of the phantom energy density, and 
must not be confused later on with the subscript with- 
out parentheses which stands for the present value of a 
quantity. 

Equation P?|) can be easily solved analytically in the 
case of power-law potentials. For a = —2 we acquire: 



(0) 



Cl e V ^"^0 - 



no / c 



C2 a 



-3/2 



(18) 



where Ci and C2 are constants determined by initial con- 
ditions. For a 7^ —2 and a (which is the trivial case 
of a constant potential), the general solution of (flT)) is 



0(0) = C'{a)a 



3/{2+a) 



(19) 



where the function C(a) is related to the potential pa- 
rameters through 



C(a) 



2a(2 -f- a)'^K 



3pmo(4 + a) 



(20) 



In (|19p we have kept only the solution part that remains 
small (together with its derivative) for small a's, in order 
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to be consistent with the matter-dominated approxima- 
tion \pff,\ <C Pm)- That is why initial-condition de- 
pendent constants are absent (which is just the central 
idea of tracker solutions). In addition, it is easy to see 
that this solution remains regular for a — > — 2. 

The zeroth-order solution for can be calculated from 
under x <C 1, that is ^^(o) = a;(o)Pm + ^(0(o)); where 



-(a0(o))^/6- The result is: 



-3p„ 



C{af 



■ a 



(21) 



'a(2 + a) 

From this expression we educe the physically interesting 
result that —2 < a < in order for /^^(o) to be posi- 
tive (additionally in this case the C(a)-definition in ((20|) 
is unambiguous). The restriction of a to negative val- 
ues, that is the consideration of normal power-law po- 
tentials instead of inverse power-law ones (contrary to 
quintessence case), was expected due to the "inverse" ki- 
netic behavior of phantoms in potential slopes. The ad- 
ditional lower bound of a values, which breaks this form 
of symmetry between phantom and quintessence models, 
arises from the extra and necessary assumption of posi- 
tivity of the total phantom energy, which is not present 
in quintessence case 0, |^, |^ where one can always ob- 
tain it (the potentials can always be shifted to positive 
values in cosmology). In other words, if the potential is 
too steep (large —a), then at early times (where the field 
value <C 1, i.e it is close to the potential minimum) the 
phantom potential energy will be very small (oc cj)^") and 
thus unable to make the total energy positive. On the 
other hand, if the potential is smooth (small —a) then 
positivity is obtained. Quantitatively, a — —2 separates 
the two regions, and indeed one can see that the energy 
density arising form is negative unless ci, C2 0, in 
which case it is exactly zero. But ci,C2 — ^ is implied 
by the requirement 0(o)i'/'(o) ^ 1 a-t small a. Therefore, 
for the limiting case a — —2, the matter-dominating con- 
dition leads to P0(o) = 0, and thus this exponent value 
bounds the two aforementioned regions. 

The corresponding zeroth-order solution for is 
P0(o) = X(_o)p,n " V^((/>(o)), which leads to: 



P4>{0) — 6pmO- 



C{a) 



■ a 2 + a 



(22) 



'a{2 + a)^ 

Therefore, for the zeroth-order solution of the equation- 
of-state parameter we obtain: 

P0(o) 2 a 

that is W(o) is constant. Note that —2 < a < leads to 
W(Q) < — 1, which is a self-consistency test of our calcula- 
tions since it is the basic feature of phantom cosmology. 
Finally, we can calculate the zeroth-order behavior of the 
phantom density parameter as fl^(^o) = P0(o)/(P0(o)+Pm), 
with Pm = prao/a^-, acquiring 



0(0) 



= 1 - 



1 - 



iC{af_ 
a{2 + a) 
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(24) 



Expressions ([191), (Ell>i (ES]) and are the tracker 
solutions for phantom cosmology in power-law poten- 
tials. Equivalently they can be expressed as a function 
of time, considering a oc t^^^ since we are in the matter- 
dominated era, as: 



oc i2/(2+a) 

,-2o/(2+a) 



(25) 
(26) 



with the equation-of-state parameter given by 
These expressions provide an excellent approximation to 
the behavior of the phantom field as long as p^, \p^\ <ti 

Pm- 



III. FIRST-ORDER PERTURBATION OF 
MATTER-DOMINATED COSMOLOGICAL 
EVOLUTION 

As cosmological evolution continues, the phantom field 
increases and phantom energy density becomes non- 
negligible, although still dominated by the dark matter 
one. Thus, we expect that the various quantities will pro- 
gressively start diverging from the expressions obtained 
above. In this section we are interested in studying an- 
alytically this intermediate evolution stage, and thus we 
perform a first-order perturbation to the zeroth order so- 
lutions of section |TT1 The corresponding quantities are 
denoting by the subscript 1 in parentheses, and the total 
ones by tilde, i.e: 

0(o)+</'(i), P0 = P0(o)+P0(i), w = W(o)+u;(i). (27) 

This perturbation is equivalent with keeping terms up to 
first order in x (i.e keeping only a;(o)) and in 0(i)/0(o), 
in the exact evolution equation (jl6p . In addition, one 

can also expand the potential as V{(t>) — V"((/)(o)) -I- 
+ C'(</'(i))- Therefore, inserting the expan- 
sions (f27l) in (fTH]) . we obtain 



2 J," 

a 0(1) 



50(1) - 3a;(o)0(o) 

^ 3p^(o) dV 
(Pmf deb 



3a0(p)y 

2pm 



(28) 



where 



[;(o) — — (a^'^Q^^/G. Thus, for the power-law po- 



tential dH) we acquire 



1 



-2- 
(0)_ 

PmO 



3a(l + a)/va30-2-" 



27Ka30i-" 
2(2 + a)p,nO 



0. 



The general solution of (|29p. for — 2 < a < 0, is 

3(6 + a) ^3 



0(1) 



a{2 + a){a^ + 8a + 28) 



(0)' 



(29) 



(30) 
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where we have kept only the part that remains small 
(together with its derivative) for small a's, in order to be 
consistent with the matter-dominated approximation. 

The perturbation of can be calculated from (fT^ 
keeping the corresponding terms, thus: 

0(1) 

P0(i) = 2^(0)^0(0) + X(i)Pm - a- — V, (31) 



where = ~3'^^'/'(o)'?^(i) ''^^ arises from x — a;(o) + 
a;(i) = — a^((/)'^Q^ + 20'^pj0|^j)/6. Therefore, using also 
(fT9|) . ([2T|) and ([30]) in order to express the result in terms 
of 0(0) and 0(i), we acquire 



P0(i) 



q(4 + g ) ^0(1) 
6 + a 



P0(O)- 



(32) 



Similarly, the perturbation of w can arise form (|15p as 

P0(1) P0(O) 



«^(i) 



P0(O) 



pI{o) 



P0( 



1) = (l-W(o))( + 

^P0(O) 0(0) ^ 

(33) 

since p^(i) = 2:(o)P0(o) + a;(i)Pm + ^J^^V. Therefore, the 
final result writes 

q;(4 + a) /0(i) \ 

Expressions ([5^ and ([55]) coincide with the correspond- 
ing one for quintessence in inverse power-law potentials 
0, however and p^(o) has a sign difference and in 
addition a is negative and bounded. Thus, 0(i) is nega- 
tive and, contrary to quintessence, p^(i) is negative while 
is positive. 

In conclusion, the total solutions for intermediate 
times, that is at low redshifts, read: 



1 



3{6 + a)C{af 



a{2 + a)(a'^ + 8a + 28) 
P4>io) + P0(i) = 



-,6/(2+a) 



P0(O) 



1 



3{4: + a)Cia) 



{2 + a)ia^ +8a + 28) 



-,6/(2+a) 



,(35) 



,(36) 



and 



w 



w 



(0) 



"(1) 



3(4 + a)C(a)2 



.6/(2+a) 



(2 + a)(a2 + 8a + 28)" J'^^^^ 

where C(a) is given by ([^0)1 . Finally, 17^ can be easily 
calculated as ilcf, = P0/(p0 + Pm), giving 

3(4 + a)C(Q)2 



= 1-n + 



P4>{0) 



(2 + Q)(a2 + 8a + 28) 



3Q + 12 
a 2 + 



(38) 

Lastly, note that since we have these quantities as a func- 
tion of the scale factor, it is straightforward to express 
them as a function of the redshift z through — = 1 + z, 
with flQ = 1 the present value. 



IV. COMPARING ANALYTICAL AND 
NUMERICAL RESULTS 

In the previous two sections we extracted analytical 
expressions for phantom evolution in power-law poten- 
tials, when the universe is dominated by the dark matter 
sector, and we argued that these solutions are valid at 
early and intermediate times, that is at both high and 
low redshift. In order to test the precision of our formu- 
lae, we evolve numerically the exact cosmological system 
calculating the exact quantities 0(a), f20(a), w{a), and 
then we examine their divergence from the zeroth and 
first order expressions by investigating the correspond- 
ing ratios. We choose initial conditions consistent with 
initial matter domination (if this requirement is fulfilled 
then the results do not depend on the specific initial con- 
ditions), and we fix k in order to acquire flmo ~ 0.28 and 



00 



0.72 at present. 



In fig. [T] we present the ratios of the zeroth (tracker) 
and first order field solutions to the exact numerical 
value, that is 0(o)/0 and 0/0, as a function of the red- 
shift. The calculations have been performed for three po- 
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FIG. 1: (Color Online) The ratios <f>(Q)/(l> (upper three curves) 
and (j)/4> (lower three curves) as a function of the redshift, for 
a — —0.5 (black, solid), a = — 1 (red, dashed), and a 
(green, dotted). 



1.5 



tential cases, namely a = —0.5, a — —1 and a =^ —1.5. 
As we observe, 0(o)/0 is very close to 1 for z > 1.5, that 
is the tracker solution is a very good approximation at 
this early evolution stage. For 1.5 > z > 0.5 the zeroth 
order solution is not a good approximation (with error 
8%) but the first order one remains within 98% accuracy 
(with even better results for larger |a| values). How- 
ever, after that stage, the phantom cosmological effects 
become significant and our approximation breaks down 
rapidly. 

Similarly, in fig. [5] we present 17^(0)7^0 and Q^/ft^ 
versus z. The divergence of the zeroth order solution from 
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FIG. 2: (Color Online) T/ie ratios r2^(o)/f20 (upper three 
curves) and 0.^/0.,), (lower three curves) as a function of the 
redshift, for a = —0.5 (black, solid), q = — 1 (red, dashed), 
and a = —1.5 (green, dotted). 
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FIG. 3: (Color Online) The ratios iu(o) /w (upper three curves) 
andw/w (lower three curves) as a function of the redshift, for 
a — —0.5 (black, solid), a = —1 (red, dashed), and a = —1.5 
(green, dotted). 



the exact one starts at z « 1.5. However, as we observe, 
the first order solution is very accurate (with error less 
than 2%) up to z « 0.3. After that point, the first order 
solution starts diverging rapidly from the exact one and 
our approximation is not vaUd. The fact that at z = 
in some cases (large \a\) the zeroth order solution seems 
to be closer to the exact evolution comparing to the first 
order one, is a result of the signs of the corresponding 
terms (namely positive (/)(o) f-i^d negative 4>(i)), but it 
has no meaning since our approximations are not valid 
in that regime. 

Finally, in fig. [3] we depict w^q^/w and w/w as a func- 
tion of the redshift. As we see, the tracker solution is very 
accurate up to z ~ 1, while the first order one agrees with 
the exact evolution within 2% up to z sa 0.3. 

In summary, we observe that the tracker (zeroth order) 
solution is accurate within an error of 2% at early times 
and up to z « 1.5. At intermediate times we have to use 
the first order analytical solution, which agrees with the 
exact cosmological evolution within an error of 2% up 
to z « 0.5. After that point, the phantom dark energy 
sector enhances, it dominates the cosmological evolution, 
and our approximation breaks down rapidly. 



V. DISCUSSION AND CONCLUSIONS 

Since we have tested the accuracy of our analytical ex- 
pressions and we have determined their applicability re- 
gion, one can use them to describe an arbitrary phantom 
evolution in power-law potentials up to z = 0.5. Here 
we will present some additional applications concerning 
observable quantities. 

Having extracted expressions for w{a) and 17^ (a) which 



are valid up to low redshifts, we can derive the corre- 
sponding w{Vl(f,) relation . Indeed, using pT|) . ([37]) and 
([55)1 we result to 

yjin.) = — — - '^-^^^^ -n.. (39) 

^ ^' 2 + a (2 + a)(a2 -K8a-K28) ^ ' 

Thus, at early times, where phantom energy density 
(that is fi^) is negligible, w is equal to a constant value, 
namely W(q) . At intermediate times it increases linearly 
with f20. Comparing with the corresponding result for 
quintessence in inverse power-law potentials [9|] there are 
two differences. Firstly, w{^l^) is an increasing function 
instead of a decreasing one. Secondly, since a lies now 
in the interval (—2,0), the slope is not a slowly-varying 
function of a but it ranges from 0.1 at a = —0.5 to 0.8 
at a = —1.5. Finally, one can examine the accuracy of 
expression p9l) by comparing it to the exact evolution. 
Indeed, we find that it is very satisfactory (with an error 
of less than 2%) up to z « 0.5 (or equivalently up to 
il^ ~ 0.3), and then the error grows reaching to w 9% at 
z = 0. 

Additionally, and similarly to 0] , we can use expression 
P7p for {y(a) in order to acquire a parametrization for 
w{z) for redshifts up to z « 0.5, where our approximation 
is still valid. Indeed we find the same result with 
namely: 

wfit = W(^o) + {wo - «;(o))(l + z)-6/(2+a)^ (40) 

with the obvious difference that now the undetermined 
parameter wq, which is just the present value of w, is 
smaller than —1. This fit for phantom evolution in power- 
law potentials proves to be satisfactory within an error 
of 3%. Finally, one can examine (HOI) in parallel with the 
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phenomenological parametrization = wq+wi ln(l + z) 
of tl2j, where he finds that the two expressions have a 
comparable level of accuracy [§| . 

In this work we have investigated phantom models with 
power-law potentials. We have extracted the tracker so- 
lutions under the assumption of matter domination, that 
is the general solution and common behavior of all such 
models at early times, that is at high redshifts. Due to 
the nature of phantom fields, which present an "inverse" 
kinetic behavior in potential slopes, we find that total en- 
ergy positivity requires normal power-law potentials in- 
stead of inverse power-law ones (contrary to quintessence 
case) , with the potential exponent being bounded by the 
quadratic form. 

Furthermore, we have extracted the general cosmologi- 
cal solution at intermediate times, that is at low redshifts, 
which is the period during the transition from matter to 
dark-energy domination. Such a solution can be very 
useful in dark energy observations, since these (super- 
novae la, WMAP and SDSS ones) are related to this 
cosmological period [l| . The comparison of the analytic 
expressions with the exact evolution arising from numer- 
ical elaboration, shows that the tracker solution agrees 
with the later within 2% for redshifts z > 1.5, while the 
first order solution is accurate within 2% up to z w 0.5. 

In the aforementioned analysis we have not used ob- 
servations in order to further restrict the potential form. 



desiring to remain general. However, from the w{a) ex- 
pression derived in (j37p we conclude that smaller values 
of I a I are favored if we want to acquire aw at z — 0.5 
which is not too negative, consistently with observations 
within 2(7. Note also that smaller |a|-values facilitate and 
enhance the total energy positivity. 

Finally, we have to make a comment about the quan- 
tum behavior of the examined model. As it is known, 
the discussion about the construction of quantum field 
theory of phantoms is still open in the literature. For 
instance in l^j the authors reveal the causality and 
stability problems and the possible spontaneous break- 
down of the vacuum into phantoms and conventional 
particles. On the other hand, there have also been 
serious attempts in overcoming these difficulties and 
construct a phantom theory consistent with the basic 
requirements of quantum field theory with the 

phantom fields arising as an effective description. The 
present analysis is just a first approach on phantom 
fields in power-law potentials. Definitely, the subject of 
quantization of such models is open and needs further 
investigation. 
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